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ABSTRACT 

Based on long baseline (5-7 years) multi-epoch HST/ ACS photometry, used previously to measure the proper 
motion of M31, we present the proper motions (PMs) of 13 main-sequence Milky Way halo stars. The sample 
lies at an average distance of r ~ 24 kpc from the Galactic center, with a root-mean-square spread of 6 kpc. At 
this distance, the median PM accuracy is 5 km s^^. We devise a maximum likelihood routine to determine the 
tangential velocity ellipsoid of the stellar halo. The velocity second moments in the directions of the Galactic 
{I, b) system are 



w,-;^/^ = 123l23 and (w?)^/^ = 83t?R km s"^ We combine these results with the 
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known line-of-sight second moment, (vf^^) ^1"^ ~ 105 ± 5 km s ^, at this (r) to study the velocity anisotropy of 
the halo. We find approximate isotropy between the radial and tangential velocity distributions, with anisotropy 
parameter /? = O.OIq ^. Our results suggest that the stellar halo velocity anisotropy out to r ^ 30 kpc is less 
radially biased than solar neighborhood measurements. This is opposite to what is expected from violent 
relaxation, and may indicate the presence of a shell-type structure at r ~ 24 kpc. With additional multi-epoch 
HST data, the method presented here has the ability to measure the transverse kinematics of the halo for more 
stars, and to larger distances. This can yield new improved constraints on the stellar halo formation mechanism, 
and the mass of the Milky Way. 



1. INTRODUCTION 

The oldest, and most metal-poor stars in our Galaxy reside 
in the stellar halo; a diffuse envelope of stars extending out 
to r ~ 100 kpc. The orbital timescales of these halo stars 
are very long compared to the age of the Galaxy, thus the 
phase-space structure of the stellar halo is intimately linked to 
its accretion history. Furthermore, the extreme radial extent 
of halo stars, well-beyond the baryonic center of our Galaxy, 
makes them excellent tracers of the dark matter halo. 

Stars diffuse more quickly in configuration space as op- 
posed to angular momentum space, so often the velocity 
structure of stellar halo stars provides the strongest link to 
their initial conditions. Global kinematic properties, such 
as the relative pressure between tangential and radial veloc- 
ity components, otherwise known as the velocity anisotropy, 
can provide important insight into the formation of the stellar 
halo. For example, most models of vi olent relaxation (e.g. 
iDiemand et al.l 120041 : ISales et akl 120071) . predict an increas- 
ingly radially biased velocity ellipsoid with distance. Local 
studies, limited to heliocentric distances Z) < 10 kpc, have 
utilized full three-dimensional (3D) kinematics of halo stars 
and find a strongly radially bias ed velocity anisotropy with 
/3 = l -fa?)/2(z;g ^ 



200llKepleY et al]|2007llSmith et al.ll2009HBond et al]l2oTo[ 



iially bias ea velocity anisotropy witn 
0. 5-0.7 (iChiba & Beersl2000l: [Gouldl 



Carolloetal.ll2oTol) . 

Beyond £) > 10 kpc, we are limited to one velocity com- 
ponent: the line-of-sight (LOS) velocity. At large distances, 
where D ^ i?o, the LOS velocity is almost identical to the ra- 
dial velocity component, and hence we have a very poor han- 
dle on the tangential motion of halo stars. Despite this short- 
coming, several studies have used large samples of halo stars, 

' Department of Astronomy and Astrophysics, University of California 
Santa Craz, Santa Cruz, CA 95064, USA; alis@ucolick;.org 
^ Hubble Fellow 

^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, 
MD 21218, USA 



widely distributed over the sky, t o tease out t he hal o velocity 
ellipsoid beyond L» 10 kpc. ISirko et al.1 (12004 modeled 
the LOS velocities of iV 1000 blue horizontal branch stars 
(BHB) selected from the Sloan Digital Sky Survey (SDSS) 
data release 4 (DR4) ranging from 10 < D/kpc < 30, and 
found a velocity ellipsoid with /? ~ —0.1 ± 0.2, which is 
consistent with isotropy. More recentlv. lDeason et al.l(l20I2ah 
used iV - 2000 BHB stars selected from SDSS DR8 be- 
tween 20 < D/kpc < 40, to simultaneously derive the ve- 
locity anisotropy and mass profile of the Galaxy: they found 
a rad ially biased veloc ity anisotropy with /? ^ 0.5. Fi- 
nall v.lKafle et alJ (l20I2h analyzed the SDSS DR8 BHB sam- 
ple jXue et al.l (1201 Ih )? and found, by dividing the sample 
into ~ 15 radial bins, that the anisotropy profile shows a 
sharp decline at r '--^ 17 kpc to /3 ~ —1 and is radially bi- 
ased (with ji ~ 0.5) either side of this apparent dip (see also 
ISamurovic & Lalovicll20I H) . 

These observations point to a fairly complex velocity 
anisotropy profile at large distances, which is likely affected 
by substructure in the stellar halo. Unfortunately, analyses 
using only LOS velocities require modeling assumptions re- 
garding the underlying potential, and it is not obvious how 
such systematics may bias the results. Clearly, in the ideal 
case, one would like to directly measure the tangential mo- 
tions of halo stars. This is a daunting task at large distances 

in the halo: at Z? ^ 10 100 kpc, a tangential velocity of 

Vt ^ 100 km s^^ corresponds to a proper motion (PM) on the 

sky of ^ ^ 2 0.2 mas yr^^. This requires an astrometric 

accuracy that is un-feasible for current stellar PM surveys. 

The Hubble Space Telescope (HST) has unparalleled astro- 
metric capabilities: the displacement of stars over time rela- 
tive to a stationary background source can be used to produce 
very accurate PMs. Quasi-stellar objects (QSOs) have con- 
ventionally been used as reference sources in HST PM stud- 
ies (e.g. iKallivay alil et al. 2006a: Kallivayalil et_aL 2006bl: 
iPiatek et al.ll2008h . In addition, distant background galaxies 
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Figure 1. Color-m agnitude diagrams (CMDs) of the three HST/ ACS M31 fields (Spheroid, Disk and Stream), using the photometry from 
iBrown et al.l ( 120091) . The black box indicates our selection of candidate foreground Milky Way halo stars. The symbols indicate the stars in 
each FOV which fall into our CMD selection box (11 in Spheroid, 9 in Disk and 3 in Stream). Blue squares are Milky Way halo stars, red 
triangles are M31 stars, and a possible Milky Way disk star is highlighted with a magenta circle (see Section l23t . 

can also be used as stationary reference sources. The advan- 
tage of background galaxies is that there are many more of 
them over the HST field-of-view, therefore a \/N averaging 
advantage means that they can potentially yield more accu- 
rate PM measurements. Galaxy positi ons can be difficult to 
measu re (compared with QSOs), but iMahmud & AndersonI 
(|2008|) have recently developed a template-fitting technique 
to measure accurate positions of ba ckground galaxies. This 
technique has been extended by (ISohn et al.l l2012allbh to 
measure accurate PMs of local group galaxies using multi- 
epoch HST photometry. This work has yielded unprece- 
dented PM accuracy (~ 0.01 mas yr^^) and has been used 
to derive the tangential motion of M31, located at a dis- 
tance of ^ 770 kpc, and by ex tension, the total mass of 
the local group (S ohn et al. 20 12atlvan der Marel et al.l2012at 
Ivan der Marel et al..,2012b.) 

In this stu dy, we make use of the unparalleled PM accuracy 
achieved bv lSohn et al.l ( l2012ah . to extract the PMs of individ- 
ual halo stars in the foreground of these HST/ ACS fields. The 
HST fields are small, but deep, which allows us to study fore- 
ground main-sequence (MS) halo stars. This is in contrast to 
spectroscopic studies, which cannot go deep enough to study 
this dominant stellar halo population. For the first time, we 
are able to directly measure the tangential motions of distant, 
field halo stars. 

The paper is arranged as follows. In Section 2 we describe 
the multi-epoch HST/ ACS imaging, and the techniques used 
to extract PMs of halo stars in these fields. In Section 3 we de- 
scribe our analysis of the halo stars, and outline our maximum 
likelihood routine used to determine the halo velocity ellip- 
soid. Section 4 describes our results, and in Section 5 we dis- 
cuss the implications of our derived halo velocity anisotropy. 
Finally, we summarize our main findings in Section 6. 



scribed inlBrown et a l.l(l2006|). and point source catalogs were 
distributed in Brown et al.l (120091) . Program GO-11684 re- 
observed these fields using HST ACSAVFC and WFC3/UVIS 
in order to measure the proper motions of their stars and M3 1 
itself. These data provides five to seven year baselines with 
respect to the pre-existing first-epoch data. The second-epoch 
observations, data red uction, and determi nation of proper mo- 
tions are described in lSohn et al.l (l20I2al) (and see below). 

2.1. Proper Motions 

During the course ofthe lSohn et"aI1(l20I2al) M31 PM study, 
PM catalogs for individual stars in the three HST fields were 
created. As the M3 1 star s were used to align first- and second- 
epoch data in the Sohn et al.l ( l2012ah study, these PMs are 
measured relative to the M31 stars. To obtain absolute PMs, 
we added the mean absolute PMs of the M31 sta rs for each 
field. For this we used the values reported by iSohn et aTl 
( I2012al) . as determined with respect to stationary background 
galaxies. The PM errors were calculated based on the posi- 
tional repeatability seen in the multiple first-epoch ACSAVFC 
measurements. We adopt similar errors for WFC3/UVIS as 
for ACS/WFC, which may be a slight overestimate since the 
WFC3/UVIS pixel scale is smaller. 

In the following section we select potential foreground 
Milky Way halo stars in the three HST /ACS fields. The proper 
motions of these sta r s are th en extracted by cross-matching 
with the ISohn et al.l (l2012al) catalog of individual star PM 
measurements. In the magnitude range under consideration 
(21.5 < mF8i4W ^ 25.5; see below), the average uncertainty 
in the proper motion measurements is (t^ ~ 0.05 mas yr^^. 
The corresponding velocity error depends on the distance, but 
for example is only 5 km s^^ at Z? = 20 kpc. This is consid- 
erably more accurate than previous proper motion measure- 



2. //ST DATA 

The data employed in our analysis come from three HST 
observing programs: GO-9453, GO-10265 (PI: T.Brown), 
and GO-11684 (PI: R.P van der Marel). Programs GO- 
9453 and GO-10265 obtained deep optical imaging of three 
fields in M3 1 (Spheroid, Disk and Stream: see Figure [T} in 
two filters (F606W and F814W) using HST ACS/WFC, with 
the primary goal of measuring the star formation history in 
these fields. The observations and data reduction are de- 



ments of individual stars in the halo (with tT„ 
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-4 mas 



yr~^; see e.g . 'Munn et alj|2004t iCasetti-Dinescu et al.ll2006l: 
iBramich et a l. 2008) 

2.2. Selection of Foreground Milky Way Halo Stars 

We select candidate foreground Milky Way halo stars us- 
ing the color-magnitude diagrams (CMDs) of the HST/ACS 
M31 fields. The main-sequence turn-off (MSTO) for Milky 
Way halo stars is located at brighter (apparent) magnitudes 
than the M31 MSTO, but at similar colors blue-ward of 
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Table 1 

The properties of candidate foreground halo stars selected from the M3 1 CMDs. We give the right ascension (RA) and declination (DEC), 
HST/ACS S TMAG magnitudes, Galactic PMs and designated class (Halo, MW Disk or M31 star). The RA, DEC and magnitudes come from 
iBrown et al . (2009), and the prop er motions derive from the study by Sohn et al. (2012a). Variables stars in M31 identifi ed b y Brown et al.j 
120041 (B04) and lJefferv et alj20lil (Jll) are indicated. One of these stars, VI 16, was also identified as an M31 member bv lKalirai et all ilOOB) 
as part of the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. 



"T-Feoew— ?TiF8i4W ^ —0.3. InFigurefT] our color-magnitude 
selection box is shown on the CMDs of the three M31 fields. 
We target halo stars in a sparsely populated region of the 
CMD, where we expect the least contamination from Milky 
Way disk stars and M31 stars. We select 13, 9 and 3 stars from 
the Spheroid, Disk and Stream fields respectively, giving a to- 
tal of 23 candidate halo stars from the three M31 fields. The 
extracted PMs and HST/ACS STMAG magnitudes are given 
in TablelU 

2.3. Comparison with Besangon Galaxy model 

In Figure [2] we compare our C MD halo star se lection with 
the Besangon Galaxy model (R obin et al.l 120031) . We select 
stars from this model in a 1 deg^ field-of-view (FOV) centered 
on M31. We convert Johnson-Cousins UBVRI photometry to 
HST/A CS STMAG using the relations given bv lSirianni et al.l 
(120051) . The stars inside our color-magnitude selection box are 
predominantly halo MS stars, and a small fraction (7 %) are 
disk stars. Most of these disk stars are intrinsically faint white 
dwarfs, which have similar (blue) colors to MSTO stars. Note 
that if we scale this 1 deg^ FOV to the HST/ACS FOV then the 
Besangon model predicts N ^ 5 halo stars per HST pointing. 
Our sample of iV = 23 stars from three pointings is slightly 
higher than this prediction; this is due to contamination by 
M31 stars in our sample, which we now discuss. 

Figure [3] shows the Galactic proper motion components of 
the Besangon Galaxy model stars which lie within our color- 
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Figure 2. A CMD for stars in the Besanjon Galaxy model. The FOV 
is 1 deg^ and is centered on M31. The black/gray shaded regions 
indicate halo stars and the magenta squares indicate disk stars (only 
25% of disk stars are shown for clarity). Within our selection box, 
7% of the foreground stars are disk stars. 

magnitude selection box. The gray crosses are halo stars and 
the magenta circles are disk stars. The disk stars have a much 
higher spread in PM than the halo stars. This is because the 
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Figure 3. Proper motions (PMs) of stars in the Besanfon model in 
Galactic coordinates {^i, iib). The stars shown cover a I deg^ FOV 
centered on M31, and fall in the CMD selection box given in Figure 
[T] The gray crosses are halo stars and the magenta circles are disk 
stars. The blue squares show the PMs of stars selected from the M31 
CMDs in Figure [T] The distribution of PMs is offset in the positive 
/i; direction; this is due to the reflex motion of the Sun. There is a 
concentratio n of (9) stars with z ero PM, consistent with the net PM 
of M31 (see ISohn et ani2012d ; Ivan der Marel et al.ll2012al) . These 
stars are highlighted with red triangles, and the dashed-yellow circle 
shows the region in PM space consistent with M31 stars. The black 
contour encompasses 99% of the halo population in the Besangon 
model. Outside of this contour, the bias caused by including high 
PM disk stars becomes non-negligible. The bottom panel is a zoom- 
in of the top panel. 

stars belonging to the Milky Way disk are much closer than 
the halo stars, so they tend to have much larger proper mo- 
tions. The black contour indicates the PM region contain- 
ing 99% of the halo stars. Within this region there is a very 
small amount of disk contamination (~ 1%), and the low PM 
disk stars have a negligible effect on our an alysis described in 
Section [3] (see discussion in Section l4.2.3b . The blue squares 
(and red triangles) show the PMs of the 23 stars selected from 
the HST/ ACS images. Encouragingly, the distribution of PMs 
closely resembles the halo population of the Besangon model. 

The dashed yellow circle indicates the region in PM space 
consistent with the net PM of M31 ((mz,m31, Mfc,M3i) = 



(0.04, -0.03) ± 0.01 mas yi'^; iSohn etaUlMUh . Note 
that this region in PM space has a radius of ~ 0.5 mas yr~^. 
There are TV = 9 stars in our sample that lie within this re- 
gion, which are highlighted by red triangles. Given the un- 
certainties in our proper motion measurements and the inter- 
nal velo city dispersion in the M31 fields {Vint ^ 125 km 
s~^, see Ivan der Marel et al.ll2012ab . we find that these nine 
stars are typically within 3-4 a of the net motion of M3 1 . In 
the following analysis we excise these stars from our sample 
as they are, presumably, horizontal branch, asymptotic giant 
branch and/or variable M31 stars. In fact, four of these stars 
have been id entified as variable stars in M3 1 from indepen- 
dent studies dBrown et al.ll2004l: iJefferv et al.ll201 ll) . 

Of the remaining iV = 14 stars, one lies outside the 99% 
contour of the Besangon model, this star is highlighted with 
a magenta circle. Given the prediction of 7% disk contami- 
nation from the Besangon model, it's likely that we may have 
^ 1 disk star in our sample of 14 stars, and the location of 
this star in PM space makes it the most probable candidate. 
In our analysis below, we only consider stars within the 99% 
contour defined by the Besangon model, and so we do not in- 
clude this star. However, we note that these two biases in PM 
space (due to M31 stars and disk stars) sue fully accounted for 
in our analysis. Furthermore, we verify that a slightly differ- 
ent choice of contour, which includes this possible disk star, 
does not change the final velocity ellipsoid results by more 
than ~ 15 km s~^, which is less than our statistical uncertain- 
ties (~ 30 km s~^). In the following section we outline how 
we model the tangential velocity components of the remain- 
ing = 13 halo stars. 

3. ANALYSIS 
3.1. Distance Calibration 

The tangential motion of the halo stars depends on their PM 
and distance: vt = 4.74047 /i D. Here, the PM is in mas yr~^ 
and the heliocentric distance is in kpc. 

The large spread in absolute magnitude of stars near the 
MSTO means that we cannot determine an accurate distance 
to individual halo stars. Therefore, we cannot measure the 
velocity ellipsoid directly, despite the fact that we have very 
accurate PMs. However, we can determine a probability dis- 
tribution function (pdf) of halo star distances, based on a suit- 
ably chosen model for the star formation history and 3D den- 
sity distribution of the halo (both of which have been well 
constrained by existing studies). This allows a statistical de- 
termination of velocity ellipsoid, without knowing the actual 
3D velocities of individual stars. 

We calibrate the pdf of halo star absolute magnitudes via a 
two step process; First, we use suitably chosen isochrones to 
describe the distribution of absolute magnitudes as a function 
of color Second, we translate this distribution into a contin- 
uous, analytic pdf t hat c an be used in a maximum likelihood 
routine (see Section U!2T i. 



3.1.1. 



Weighted Isochrones 

We use the lVandenBerg et"an (12006 ) isochrones, calibr ated 
for HST/ ACS STMAG photometry ferown et all 120051) . to 
model the absolute magnitudes of our selected halo stars. 
In Figure |4| we show isochrones for a range of metallicities 



ISohn et al] 42012aD gave values in the (W,N) frame. However, at the 
location of M3 1 , the (W,N) and (l,b) coordinate frame are aligned to within 2 
degrees. 
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Figure 4. IVandenBerg et alj ( 120061) isochrones in the HST/ACS 
STMAG system. Metallicities in the range -2.3 < [Fe/H] < -1.0 
are shown by the colored Unes, and ages between 10 < T/Gyr < 14 
are shown by different line styles. All these isochrones have an alpha 
to iron-peak element ratio [a/Fe] = 0.3. 

(-2.3 < [Fc/H] < -1.0) and ages (10 < T/Gyr < 14) ap- 
plicable for a halo population. In the color range of our sam- 
ple, tof606W^™F814W < —0.3, the halo stars are close to the 
MSTO, but they can also lie at fainter magnitudes on the MS, 
or brighter magnitudes on the red giant branch (RGB). Thus, 
the absolute magnitudes can range from 3 < A/f814W S '''•5, 
and our sample of halo stars potentially spans a distance 
range, 10 < L>/kpc < 100. 

We model the absolute magnitudes of the halo stars as a 
function of color, and apply three weighting factors to the 
isochrones: 

• Initial-mass-function (IMF) weight: We assume a 
Salpeter IMF and weight the contribution from differ- 
ent mass stars accordingly. Low mass MS stars thus 
have higher weight than higher mass RGB stars. 

• Metallicity weight: We assume the halo population 
has a Gaussian distribution of metallicities with mean 
[Fe/H] = —1 .9 and dispersion a = 0.5 (see e.g. 
IXue etal.1 '2008'). We assume the metallicity distribu- 
tion is constant with radius; this is in good agreement 
with recent studies which find no significant metallicity 
gradient in the stellar halo (Ma Z. et al. in preparation). 

• Age weight: We weight the different age isochrones by 
assuming a Gaussian distribution of ages with mean 
(r) = 12 G yr and dispersion <t = 2 Gyr (see e.g. 
lKaliraill20T2l) . 

In Fig. |5]we show the absolute magnitude distribution in 
six color bins, weighted by IMF, metallicity and age. The 
dashed blue line shows a double-Gaussian fit to each color 
bin. At redder colors the distribution of absolute magnitudes 
separates into the MS and RGB populations, while at bluer 
colors the populations merge into the MSTO. 

3.1.2. Analytic absolute magnitude calibration 

In Figure |6] we show the parameters from the double- 
Gaussian fits (amplitude, mean and sigma) for both the RGB 
and MS populations as a function of color. The dashed red 



lines are polynomial (second or third order) fits to the rela- 
tions. We use these polynomial relations to define a contin- 
uous, double-Gaussian pdf for absolute magnitude as a func- 
tion of color: 

G(M'f814w|™F606W— ™F814w) = Gi (Ai , Mi , CTi , Mf814w) 

+G2(A2,M2,cr2,MF814w) (1) 



(x-M)V(2a2) 



and 



where, G{A, M, a,x) = A cxp 

A, M and a (amplitude, mean and sigma) are polynomial 
functions of tof606W — '71f814W color. 

The red squares in Figure 5. show the absolute magnitude 
distributions reproduced by our analytic model. These are in 
good agreement with the (non-analytic) distributions derived 
from the weighted isochrones. 

In the following section, our absolute magnitude pdf is used 
in a maximum likelihood routine to determine the tangential 
components of the halo velocity ellipsoid. 

3.2. Maximum Likelihood Method 

In this section, we outline the maximum likelihood method 
used to derive the tangential velocity ellipsoid components of 
the Milky Way halo. We assume that the velocity distribution 
in the plane of the sky follows a Gaussian distribution, with 
constant (cri,ab,vi^o,Vbfi) over the radial range spanned by 
our data: 



Fv{vi,Vb) cx exp 



2af 



cxp 



jvb - Vb.oY 
2al 



(2) 

In the halo, where _Ro — 8.5 kpc, the Galactic co- 

ordinates, I, b are good approximations to the spherical coor- 
dinates, (p, 9 respectively. For each star, there are six observ- 
ables: the angular position on the sky (/, 6), the PM {fii, fii,), 
and the two photometric magnitudes (tof814W, ™f606w)- 
Observed heliocentric velocities are converted to Galactocen- 
tric ones by assuming a circular speed of 240 km s~^ (e.g. 
iReid et al.ll2009l:lMcMillanll20Tll:lSchonrichll2012l) at the po- 
sition of the sun (Rq =8.5 kpc) with a solar peculiar motio n 
(U, V, W)=(ll.l, 12.24, 7.25) km s^^ (!Sc honrich et al.l2010l) . 
Here, U is directed toward the Galactic center, V is positive 
in the direction of Galactic rotation and W is positive towards 
the North Galactic Pole. In the direction of M31, the velocity 
of the sun projects to: = (—139.5, 83.7) km s^^. 

The halo pdf at fixed mpeoew — ?7^f814W color, in incre- 
ments of absolute magnitude, apparent magnitude. Galac- 
tic PM and solid angle, F{y), where y is defined as y = 
i/(il/F8i4W, n^F8i4W, Mi: Mb- ^^). is given by: 



F Ay (xFypD'^G Ay 



(3) 



Here, D ~ _D(A/f814W, "^F814w) is the heliocentric dis- 
tance, Fy ~ Fy{D, ni^b) is the velocity distribution func- 
tion given in Equation|2l p = p{D, b) is the density distri- 
bution of halo stars, G = G(A/f814w|™f606W — "t-F814w) 
is the absolute magnitude pdf defined in Equation [T] and 
Ay = AA/f814wAtof814wA/j,; A^fcAfJ is the volume ele- 
ment. We assume the hal o stars follow t he brok en power-law 
density profile derived by iDeason et al.l (1201 Ibl) . but we com- 
ment on the effect of the de nsity profile parameterization on 
our results in Section 14.2.21 We marginalize over the abso- 
lute magnitude coordinate, F = /FdAfF8i4w, and define 
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Figure 5. The Afpgi4w absolute magnitude distribution of the lVandenBerg et al.l ( 120061) isochrones in six color bins. The MSTO dominates at 
bluer colors, while at redder colors the MS and RGB populations separate into two distinct populations. The blue lines show double-Gaussian 
fits to the distributions. The red squares indicate the distributions reproduced by our analytic model (see Equation[T) in these six color bins. 




Figure 6. The Gaussian parameters (amplitude, mean and sigma), which define the RGB and MS absolute magnitude calibrations, are shown as 
a function of color. The dashed-red lines show continuous polynomial fits (second or third order) to the relations. 



7 




moment^ in the plane of the sky are: (Tb 



100 -50 
V,,o [km/s] 



Figure 7. Likelihood contours for the halo velocity ellipsoid in the 
plane of the sky. The shaded gray regions and solid lines show the 
Icr and 2(7 confidence intervals, and the black crosses indicate the 
maximum likelihood parameters. 



the HkeHhood function: 



L = 



]^F(cr;,(Tf,,T;;^0,'yf),0,x) 



(4) 



Here, x denotes the observables (mF8i4w, ?tif606W ^ 
'7T.F814W, M/, Mh, b), and ai^ai,, vi fi, vt.o are free parameters. 
We use a brute force grid method to find the maximum like- 
lihood parameters. The net motion in Galactic latitude, W(,.o- 
is set to zero, while we allow for net motion in Galactic lon- 
gitude, w;,o ^ —Vci)fi, which approximates the net rotational 
velocity of the halo. 

The kinematical biases that we have introduced (see Figure 
[3]) are accounted for through a renormalization of the likeli- 
hood function over the region in observable-space that passes 
our cuts. In Equation ID F{p, x) is the normalized probabil- 
ity that a star exists at observables x, given model param- 
eters p = {ai,ab,vifi,Vbfi). We define i^sci(x) as our se- 
lection function, which states the probability that a star with 
observables x is in our sample. Thus, Fsci(x) is 1 or 0, de- 
pending on our adopted cuts. The normalized probability 
that a star is found in our sample, at observables x, given 
model parameters p is given by: F{p, x)Fsoi(x)/Af(p). Here, 
N{p) = J F(p, x)Fsoi(x) is the normalization factor, and Xc 
denotes the region in x-space that passes our selection cuts. 
We calculate these integrals numerically. This procedure re- 
duces the likelihood to: 



n 



Nip) 



(5) 



1 for any star i = 1, ...,N that 



since, by definition, Fsc\ (x; ) 
makes it into our sample. 

In the following Section, we discuss the results of apply- 
ing this maximum likelihood algorithm to our sample of halo 
stars. 

4. HALO VELOCITY ANISOTROPY 
4.1. Results 

In Figure [T] we show the results of our maximum like- 
lihood analysis (see also Table |2]i. The resulting velocity 



83 

+28 



and vifi = ^75_29 km s 



+-km 

1 



We calculate the mean and root-mean-square spread of the 
halo star distances using the pdf of our favored parameters, 
e.g. (D) = 1/N / F,niDdD, where F,„; denotes the 
pdf assuming our maximum likelihood parameters. We find 
(D) = 19 ± 1 kpc and (T]j = 6 ± 2 kpc. In Galactocentric 
radii, this corresponds to r ~ 24 kpc. 

Our maximum likelihood velocity ellipsoid quantities are 
converted into a spherical polar coordinate system using a 
Monte Carlo method. We generate stars which follow the stel- 
lar halo density profile dDeason et al. 2011b), and assign vi 
and Vb velocity distributions drawn from our likelihood con- 
straints on (T/, CTf, and vifl. In the approximate distance range 
of our sample, 10 < D/kpc < 30, the LOS velocity dis- 
tribution of stellar halo star s, is appr oximately Gaussian with 
o-ios = 105 ±5 km s-i (see lSirkoet al. 2004; Xue et al. 2008} 
iBrown et al.ll201C)l) : we apply this LOS velocity distribution 
to our Monte Carlo samples. The Vr, vo,v^ velocity compo- 
nents are then calculated from the generated star positions and 
v\os: Vb, vi velocities. In Table[2]we show the resulting veloc- 
ity ellipsoids in spherical polar coordinates. As expected, the 
vg, V(f, velocity distributions closely follow the Vb, vi coordi- 
nates. 

To calculate the velocity anisotropy we consider the ratio 
of tangential and radial pressures, i.e. {V^)/{Vr)' where 
(V^) = 0-0 + (w^)^ + cr|. The resulting velocity anisotropy 
suggests isotropic velocity pressure, with {V^)/{V^) = 
l.Olg 1 or /3 = 0.0lo;4, this is in contrast to the highly 
radial values found in the solar neighbo rhood: /3 ^ 0.7 
or iar,ad,,cre) - (140,8 0,80) km s"! ( ISmith et al.l 120091: 
iBond et alJl20I0tlCarollo e t al. 2010|). However, several stud- 
ies estimating velocity anisotropy from LOS velocity alone, 
in the distance range, 10 < D/kpc < 30, also find a 
more tangential velocity ellipsoid than the solar neighbor hood 



, .. . . _ ilgtl 

dSirko et al.l2004l;lDeason et al.l20I latlKafle et al.l20I2l) . We 
note that a significant contribution to the tangential pressure 
comes from angular momentum, as ^ 70 km s~^. In 
the following section, we discuss some potential systematics 
which could bias our results. Encouragingly, we find that sys- 
tematics play a minor role in the total error budget. 

4.2. Possible Systematics 
A.IA. Stellar Isochrones 



In Section [TT] we use the iVandenBerg et al.l (120061) 
isochrones to calibrate our halo star distances. This analysis 
relies on our assumptions of metallicity and age for the stellar 
halo. We adopt isochrones applicable to an old (T = 12 ± 2 
Gyr) and metal-poor ([Fc/H] = —1.9 ± 0.5) stellar pop- 
ulation, but it is worth remarking how these assumptions 
may affect our results. Older and/or more metal-iich pop- 
ulations have fainter absolute magnitudes and hence higher 
weight is given to smaller distances in our analysis (and there- 
fore, lower velocities). For example, the isochrones in the 
Besanfon model assume a similar halo metallicity distribu- 
tion ([Fc/H] = —1.78 ± 0.5), but assume an older 14 Gyr age 
population. If we adopt this absolute magnitude calibration 
in our likelihood analysis, then our velocity ellipsoid parame- 
ters are decreased by (Act/, Ao-f,) = (—3, —5) km s^^. Sim- 



^ We also applied our likelihood method setting [vi ) 
yields higher cr^, but the same answer for ct? + (s;^)'^ 



-(d^> =0. This 
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{vi) = -75tll 


Spherical polars 


(t;2)i/2 = 107l^ 










Velocity Anisotropy 








/(f?> _ 1 n+O-2 / - 1 A+0.4 








Position 






I = 121° 


h = -21° (£») = 19 ± 1 ± 6 kpc 


(r) = 24 ± 1 ± 6 kpc 



Table 2 

Summary of our main results. We give the velocity ellipsoid in Galactic and spherical coordinate systems and the resulting velocity anisotropy. 

Note, the LOS velocity dispersion is estimated by previous studies in the literature (see text for more details). We also give the approximate 
location of our three HST fields in the plane of the sky, as well as the average heliocentric and Galactocentric distances for our sample. For the 
latter quantities we list two uncertainties, the first being the error in the mean, and the second being the root-mean-square spread of the sample. 



ilarly, younger and/or metal-poorer populations will slightly 
increase the estimated velocity dispersions. We find that vari- 
ations in the mean metallicity and/or age by 0.5 dex and 2 
Gyr respectively, at most can change the velocity dispersions 
by ~ 10 km s~^, which is well below our observational ran- 
dom errors. We note that more substantial variations require 
significant changes in metallicity and/or age, contrary to ob- 
servational constraints and theoretical predictions. 

4.2.2. Density Profile 

We adopt a broken power-law stellar halo density pro- 
file, with halo flattening q ~ 0.59, ellipsoidal break radius 
rjj — 27 kpc , and inner and outer power-laws, = 2.3, 
ctont = 4.6 (iDeason et al.l l20lTbl) . The uncertainty in the 
halo star distances means that the density profile component 
of the pdf could potentially affect our derived velocity ellip- 
soid. For example, a shallower density profile will give higher 
weight to larger distances (and hence higher velocities), than 
a steeper density profile, and the derived velocity dispersions 
will therefore be higher. Encouragingly, recent studies show 
general agreement on the for m of the stellar halo density pro- 
file out to 50 kpc (see e.g . IWatkins et al.ll2009t ISesar et all 
120 1 U iDeason et al]|20 1 1 bh . Fo r example, if we in stead adopt 
the density profile derived by ISesar et aTl (1201 ll) (q — 0.7, 
rty = 28 kpc, Qfin = 2.6, aout = 3.8), our velocity ellipsoid 
values are increased by (Act;, Act;,) = (6, 3) km s^^. These 
are very small differences relative to the statistical uncertain- 
ties in our results (^ 30 km s~^). Note that a very steep halo 
density profile (a > 6) is required to significantly decrease 
our tangential velocity ellipsoid values. 

We also note that the systematics primarily associated with 
distances, affect our estimate of (V^^) more than the ellip- 
soid ratio (V^) / (Vg). This is because any systematic in our 
distance calibration, affects both tangential components by a 
similar degree. 

4.2.3. Disk Contamination 

The Besangon Galaxy model predicts that the majority of 
stars selected from our color-magnitude cut are halo stars, but 
~ 7% of the stars could belong to the Milky Way disk. Our 
further restriction in proper motion space (see Figure |33, re- 
duces this level of contamination even further (to 1% in the 



Besan9on model). Even so, we have quantified any bias due 
to the remaining disk stars. We select TV = 13 stars at random 
from the Besan^on model using our color-magnitude and PM 
selection, and apply our maximum likelihood analysis. We 
repeat this process using = 13 pure halo stars from the 
Besan9on model. After lO'' trials, we can quantify any biases 
due to contamination by the small number of disk stars. We 
find that any biases are negligible (< 5 km s^^), and can be 
safely ignored in the current analysis. However, we note that 
if we apply a less restrictive contour in proper motion space 
(see Figure [3]i, then these biases can increase. 

4.2.4. Solar Motion 

In recent years, the solar motion — and in particular, the 
azimuthal velocity of the Sun — has come under renewed 
scrutiny. There has been some debate over whether the cir- 
cular velocity at the position of the Sun (Rq = 8.5 kpc) is 
close to the lAU-recommended value, Vc = 220 km s~^, or 
if it needs to be revised upwards to Vg = 240 km s~^ (e.g. 
Reid et al.ll2009l:lMcMiUan|[20TltlSchonrichl2012l) . Recently, 
Bow et al.l (120 12h showed that the azimuthal velocity of the 
sun is close to ^ 240 km s~^. However, the authors claim 
that this is due to an offset between the azimuthal velocity 
of the Local Standard of Rest and the local circular velocity, 
rather than a change in the circular velocity. In either case, the 
net solar motion — the quantity of importance in this study — 
remains uncertain. 

In the context of the present study, variation of the az- 
imuthal solar motion by ^ 20 km s^^ can change the pro- 
jection of the solar motion in the direction of M3 1 by approx- 
imately (|Aw(|, Aub|) ^ (10,5) km s^^. This will most no- 
tably affect our estimate of net rotation in the halo by ~ 10 
km s^^. Thus, even with a decrease in the Sun's rotational 
velocity of 20 km s^^, we still find a significant signal of ro- 
tation in our sample. We note that the velocity anisotropy is 
largely unaffected by this translation. 

5. DISCUSSION 

Our measurement of an isotropic velocity anisotropy at 
r ^ 24 kpc, is at odds with theoretical predictions of vio- 
lent relaxation, where /3 tends to increase with radius: Our 
result oi P = O.O^q ^ is lower than solar neighborhood mea- 
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r [kpc] 



Figure 8. Velocity anisotropy as a o f funct ion radius. The co lored error bars show the estimates from nearby halo stars with full 3D veloc- 
ity measurements jBond et ani2010l: BIO; Smith et ani2009l: S09) , and distant {D > 10 kpc) halo stars with LOS velocity measurements 
dSirko et al.ll2004l : S04: lDeason et al.ll2012al : D12: lKafle et alj|2012l : K12). The red error bar shows our measure from iV = 13 halo stars with 
PM measurements in the radial range 18 < r/kp c < 30. The blue sh aded region indicates the approximate radial range of the 'break radius' 
in the Milky Way stellar halo (16 < r/kpc < 26: iDeason et ai]|2011bl) 



surements by ~ 3(7. In this section, we discuss some possible 
explanations for this surprising result. 

5.1. Global Substructure: A Galactic Shell? 

In Figure [U we show velocity anisotropy as a function of 
radius, based on various estimates from the literature. When 
our results are compared to solar neighborhood measures of (3 
at smaller radii, and LOS velocity estimates at larger radii, we 
begin to see that the velocity anisotropy may have a 'dip' at 
r ~ 20 kpc, rather than a constant or continuous decline. Fur- 
thermore, the location of this possible dip is coincident with 
a br eak in the stellar halo density: [P eason et al. (201 lb) (see 
also lSesaret"ani2011l: IWatkins et aTll2009l) find that the stel- 
lar halo density profile declines more rapidly beyond a break 
radius of r ~ 1 6 — 26 kpc. The ori gin of this break radius 
is still uncertain: iDeason et al.l (120 13h recently suggested that 
the break radius in the Milky Way may be due to a shell-type 
structure built up fr om the ag gregation of accreted stars at 
apocenter However, iBeers et al. (2012) claim that the break 
radius signifies a transition between two halo populations of 
different origin (see discussion below). In the former sce- 
nario, the kinematic signature of a shell — a predominance 
of tangential motion at the turn-around radius, plus a cold ra- 
dial velocity dispersion — is intriguingly similar to the appar- 
ent 'dip' in velocity anisotropy that we observe. At present, 
the evidence for a rise again in velocity anisotropy beyond 
r ~ 30 kpc, relies solely on LOS velocities, where the derived 
anisotropy is intimately linked to the assumptions of the halo 
potential. Therefore, it is important to independently deter- 
mine the tangential motion of halo stars beyond the break ra- 
dius in order to confirm/falsify the apparent dip in anisotropy. 



It is possible that some of our halo stars belong to the 
"TriAnd" overdensity. This structure subtends an area of 
at least 50° x 60° in the constellations of Triangulum and 
Andromeda, and is located at a si milar distance to our halo 
sample (16 5, D/kvc < 25. iRocha-Pinto et alJ 120041: 



iMajewski et ai.H2004t iMartin et al 



I2007h. The nature of 



TriAn d has bee n debated in the litera ture; iRocha-Pinto et alj 
(120041) (see also lMaiewski et aT]|2004l) propose either a bound 
core of a very dark matter dominated dwarf or a portion of a 
tidal stream, while Johnston et al. (2012) suggest that TriAnd 
is an apogalactian piece of a disrupted dwarf galaxy. The con- 
clusions we draw here are most consistent with this last inter- 
pretation. 

5.2. Stellar Halo Formation Mechanism: In-Situ Stars? 

In recent years, an additional formation mechanism for halo 
stars has been put forward. Studies based on hydrodynamical 
cosmological simulations suggests that halo stars can form 
in-situ from gas in the parent galaxy in addition to forma- 
tion in external dwarf g a laxies and then su bsequent accre- 
tion dZolotov et all 120091; iFont et al.ll20Tlb . [McCarthy et"ai1 
showed that these in-situ stars can have significant pro- 
grade rotation and therefore increased tangential pressure sup- 
port from angular momentum. Therefore, this alternative for- 
mation mechanism for halo stars could also explain the more 
tangentially biased halo star orbits. However, this does not 
explain why the orbits of halo stars in the solar neighborhood 
— presumably where in-situ stars are more dominant over ac- 
creted stars — have such strongly radial orbits. 

5.3. A Cold Stream? 



10 



Finally, our sample of halo stars is small (N = 13) and 
covers a narrow FOV, so we cannot reject the possibility that 
our results are biased by local substructure. However, we do 
not find any obvious clustering of our sample in PM space, so 
it seems unlikely that the majority of our halo stars belong to 
a common stream. Furthermore, the agreement of our results 
with studies using LOS velocities in a similar radial regime 
but different directions on the sky, suggests that the isotropic 
velocity anisotropy derived here is a true global property. 

5.4. Halo Rotation 

We find a significant signal of prograde rotation in our halo 
sample, with (v^) ^ 70 km s~^. As mentioned above, this 
signal could be due to the presence of a shell-type structure in 
this radial regime, or the halo could have a net rotation — per- 
haps due to the influence of stars born in-situ. Previous work 
on th e rotation of halo popu latio ns have foun d contrasting re- 
sults. iFrenk & Whitg (119801) and |Zinn| ([1985!) found evidence 
of net prograde rotation in the Milky Way globular cluster 

population (with (v^ ) ~ 50 60 km s^^). More recently, 

iDeason et all (12011 ah find that BHB stars with [Fe/H] > -2 
have a net streaming motion of {v^) ~ 50 km s~^ in asim- 
ilar radial range to this study. However, ICarollo et al.l (120101) 
claim that their sample of outer (D > 15 kpc) halo stars have 
a net retrograde rota tion with ^ —80 km s^^ (but see 
ISchonrichetani2011h . 

It is clear that a confused picture surrounds the rotational 
properties of the halo. This is something we hope to address 
in the future with additional HST fields. 

5.5. Future Work 

The above discussion illustrates the importance of measur- 
ing the radial velocity anisotropy profile in order to under- 
stand the formation mechanism of the stellar halo. We hope 
to extend this present study to multiple HST fields with long- 
baseline multi-epoch photometry. This will enable us to move 
beyond small number statistics, cover a wider area on the sky, 
and probe further out into the stellar halo. Studies using LOS 
velocity measurements rely on assumptions of the halo poten- 
tial and are limited to within r < 30 — 40 kpc: using PMs, we 
can more directly measure the tangential motion of halo stars. 
Furthermore, the depth of the HST fields allows us to study 
MS stars; the dominant population of the stellar halo. Current 
spectroscopic studies and the upcoming Gaia mission (with 
magnitude limit V < 20) are limited to intrinsically bright 
halo stars, such as BHB or RGB stars, which may not be un- 
biased tracers. 

Finally, an independent measure of velocity anisotropy is 
vital in order to deri v e the rn ass profile of our Galaxy. Re- 
cently, iDeason et alj ( l2012bl) found that the radial velocity 
dispersion of halo stars declines rapidly at large radii. At 
present, the degeneracy between tracer density, anisotropy 
and halo mass cannot be disentangled. However, a measure of 
the tangential motion of these distant halo stars will allow us 
to address whether or not this cold radial velocity dispersion is 
due to a shift in pressure from radial to tangential components. 
Our method, using multi-epoch HST images, may be the only 
way to measure velocity anisotropy at these large distances. 
The upcoming Gaia mission will measure proper motions for 
an unprecedented number of halo stars with V < 20, and 
will likely revolutionize our understanding of the inner stellar 
halo. However, even with bright halo tracers (e.g. BHB or 
Carbon stars), the PM accuracy of Gaia, ~ 0.3 mas yr~^ 



((TV ^ 140 km s^^ at D ^ 100 kpc), will be unable to ac- 
curately constrain the tangential motion of very distant halo 
stars. 

The PM accuracies achievable with just a few orbits of HST 
time (see Table [T]), are well below the halo velocity disper- 
sion, even at _D = 100 kpc. Therefore, PM accuracy is no 
impediment at all for studies of this kind. However, the small 
FOV of HST means that there are only a handful of halo stars 
per field. Therefore, to find the rarer stars that are at larger 
distances than the mean for our sample ((r) ^ 24 kpc), many 
fields need to be observed. 

6. CONCLUSIONS 

We derive proper motions (PMs) for iV = 13 Milky Way 
halo stars from three HST/ ACS fields with long baseline (5- 
7 years), multi-epoch photometry. The unprecedented accu- 
racy of these PM measurements (cr^ ^ 0.05 mas yr^^) allows 
us to derive the tangential velocity ellipsoid of halo stars be- 
tween 18 < r/kpc < 30. Until now, measurements of the 
tangential motion of halo stars have been restricted to within 
D < 10 kpc. This study is the first step towards measuring 
the halo velocity ellipsoid of the stellar halo beyond the solar 
neighborhood, independently of any assumptions regarding 
the underlying halo potential. We summarize our conclusions 
as follows: 

(1) We select = 23 candidate halo stars from three 
HST/ ACS fields centered on M31. These stars are selected 
from the halo MSTO region of the CMD. The individual 
PMs are extracted from the multi-epoch HST photometry. 
Inspection of the proper motions of these stars shows that 
nine belong to M31, and one is a possible Milky Way disk 
star The selected halo stars have magnitudes ranging from 
21.5 < TOF814W < 24.5 and are MS and/or RGB stars near 
the MSTO. The distances to these stars are calibrated from 
stellar isochrones applicable to an old, metal-poor halo popu- 
lation. 

(2) We devise a maximum likelihood routine to derive the 
tangential velocity ellipsoid of these halo stars. We measure 
PMs for individual stars, but not their distance or radial veloc- 
ities. Therefore, our determination of the velocity ellipsoid is 
statistical in nature. We find velocity dispersions at, = 83^ jg 
km s^^, fT; = 94j^^g km s^^ and mean streaming motion, 

vi_o ~ "75I29 '^'^ ^'^ '•^^ distance range of the halo 
star sample, 18 < j'/kpc < 30, the Galactic velocity com- 
ponents closely approximate the spherical coordinate system 

(w0 -vi and ve Vb). 

(3) The tangential velocity components are combined with 
independent measures of the radial velocity dispersion {<Tr ^ 
<T\os = 105 ± 5 km s~^) in this radial regime to estimate the 
velocity anisotropy of the stellar halo. We find approximate 
isotropy between radial and tangential velocity second mo- 
ments with l3 = O.O^Q 4- This is in contrast to the strongly 
radial anisotropy of halo stars found in the solar neighbor- 
hood {f3 ~ 0.7). The increased tangential relative to radial 
pressure has a significant contribution from angular momen- 
tum as (vcf,) ^ 70 km s^^. We note that the more tangentially 
biased velocity anisotropy outside of the solar neighborhood 
is mainly due to a decrease in radial pressure (from 140 to 
100 km s~^), rather than a significant increase in tangential 
pressure. 

(4) The radial anisotropy profile is poorly constrained, es- 
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pecially at large radii. However, there is growing evidence 
that there may be an isotropic/tangential 'dip' in velocity 
anisotropy in the radial range 15 < r/kpc < 25. Intrigu- 
ingly, this coincides with a break in the stellar density profile, 
beyond which the stellar density falls off more rapidly. These 
two lines of evidence suggest that there may be a shell-type 
structure in this radial regime. However, at present we cannot 
discount the influence of halo stars formed in-situ which may 
have more tangentially biased orbits than accreted stars. 
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